Quantum stochastic equation for the low density 
hmit 



(N 



^ ; L Accardif, A N PechenfH and I V Volovich| 

' t Centre Vito Volterra, Universita di Roma Tor Vergata 00133, Italia 

■ J Steklov Mathematical Institute of Russian Academy of Sciences, Gubkin St. 8, 

Ph ; GSP-1, 117966, Moscow, Russia 

E-mail: accardiSvolterra . mat . uniroma2 . it, pechenOmi . ras . ru, 
' volovich@ini.ras.ru 

(N ■ 

' Abstract. A new derivation of quantum stochastic differential equation for the 

^ , evolution operator in the low density limit is presented. We use the distribution 

approach and derive a new algebra for quadratic master fields in the low density limit 
by using the energy representation. We formulate the stochastic golden rule in the low 
OO i density limit case for a system coupling with Bose field via quadratic interaction. In 

' particular the vacuum expectation value of the evolution operator is computed and its 

. exponential decay is shown. 

> 

' Submitted to: Journal of Physics A: Mathematical and General 

O^' PACS numbers: 02.50.EY, 02.50.FZ, 03.65.YZ, 05.10.Gg and 42.50.Lc 

> 



§ permanent address: Steklov Mathematical Institute of Russian Academy of Sciences, Gubkin St. 8, 
GSP-1, 117966, Moscow, Russia 



Quantum stochastic equation 



2 



1. Introduction 

There are many studies of tlie large time behavior in quantum theory. One of the 
powerful methods is the stochastic limit. Many important physical models have been 
investigated by using this method (see []I], 0, ^ for more discussions). This method, 
however, is restricted to the studying of the long time behavior of the models in the 
weak coupling case, i.e. when coupling constant is a small parameter, and it can not 
be applied directly to important class of models which contain terms in the interaction 
without small coupling constant. 

This last class of models includes models in which the small parameter is the 
density. Such models naturally arise in the low density limit (LDL). The LDL for a 
classical Lorentz gas is the Boltzmann-Grad limit. For classical systems there has been 
a considerable progress in the rigorous derivation of the Boltzmann equation. Lanford 
using ideas of Grad |P proved the convergence of the hierarchy of correlation functions 
for a hard sphere gas in the Boltzmann-Grad limit for sufficiently short times. This proof 
was extended by King to positive potentials of finite range. The limiting evolution 
of the one particle distribution is governed by the non-linear Boltzmann equation. 

The test particle problem was studied by many authors(see a review of Spohn |]^). 
One considers the motion of a single particle through an environment of randomly 
placed, infinitely heavy scatterers (Lorentz gas). In the Boltzmann-Grad limit successive 
collisions become independent and the averaged over the positions of the scatterers the 
position and velocity distribution of the particle converges to the solution of the linear 
Boltzmann equation. 

To describe a quantum physical model to which LDL can be applied let us consider 
an N-level atom immersed in a free gas whose molecules can collide with the atom; the 
gas is supposed to be very dilute. Then the reduced time evolution for the atom will be 
Markovian, since the characteristic time ts for appreciable action of the surroundings 
on the atom (time between collisions) is much larger than the characteristic time in for 
relaxation of correlations in the surroundings. Rigorous results substantiating this idea 
have been obtained in [Q. 

It is known that the dynamics of the N-level atom interacting with the free gas 
converges, in the low density limit, to the solution of a quantum stochastic differential 
equation driven by quantum Poisson noise. Indeed, from a semiclassical point of view, 
collision times, being times of occurrence of rare events, will tend to become Poisson 
distributed, whereas the effect of each collision will be described by the (quantum- 
mechanical) scattering operator of the atom with one gas particle (see the description 
of the quantum Poisson process in |]T0[). 

In this paper a derivation of the quantum stochastic differential equation for the 
evolution operator in the low density limit is presented. The equation obtained is 
equivalent to the stochastic equation which has been derived in ^ but we use a new 
method. We use the distribution approach P, || and derive a new algebra for quadratic 
master fields in the low density limit by using the energy representation. An advantage 
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of this method is the simphcity of derivation of quantum stochastic equations and 
computation of correlation functions. We formulate the stochastic golden rule in the 
low density limit case for a system coupling with Bose field via quadratic interaction. 
In particular the vacuum expectation value of the evolution operator is computed and 
its exponential decay is shown. 

Main results of the paper are quantum stochastic differential equation (0), the new 
algebra of commutation relations for the master field (Theorem |^) and the expression 
for the expectation value of the evolution operator (pSl) . 

In this paper we obtain unitary evolution which is given by the solution of 
the quantum stochastic differential equation. Using this equation one can obtain 
corresponding quantum Langevin and master equations. 

An important problem in theory of open quantum systems is the rigorous derivation 
of quantum Boltzmann equation from microscopic dynamics. Aa approach to the 
derivation of classical and quantum Boltzmann equations, based on BBGKY-hierarchy, 



has been presented in the work of Bogoliubov ||rT| . The low density limit for the model 
under consideration, with completely different methods, based on quantum BBGKY 
hierarchy has been investigated by Diimcke 0. 

Let us explain our notations. We consider a quantum model of a test particle 
interacting with a reservoir (heat bath) . We shall restrict ourselves to the case of a Boson 
reservoir in this paper. Let Tis be the Hilbert space of the system (test particle) with the 
system Hamiltonian H^. The system Hilbert space for the iV- level atom Tis = C^. The 
reservoir is described by the Boson Fock space r(7^i) over the one particle Hilbert space 
Til = L^(R'^), where = 3 in physical case. Moreover, the Hamiltonian of the reservoir 
is given by iJ^ := dr(iJi) (the second quantization of the one particle Hamiltonian Hi) 
and the total Hamiltonian of the compound system is given by a self-adjoint operator 
on the total Hilbert space Tis ® r(7ii), which has the form 

Htot := Hs^l + l® Hk + Hint =■ H^,^^ + H.^f 

Here i^int is the interaction Hamiltonian between the system and the reservoir. The 
evolution operator at time t is given by: 

Obviously, it satisfies the differential equation 
dtUt = -iHUt)Ut 

where the quantity Hi^itit) will be called the evolved interaction and defined as 

The interaction Hamiltonian will be assumed to have the following form: 
Hir^t -.= 0® A+((7o)A((7i) + ® A+igi)Aigo) 

where Z) is a bounded operator in Tis, D G 6(7^^), A and are annihilation and 
creation operators and g^^gi G Tii are form-factors describing the interaction of the 
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system with the reservoir. This Hamiltonian preserves the particle number of the 
reservoir, and therefore the particles of the reservoir are only scattered and not created 
or destroyed. This model was considered by Davies [12| in the analysis of the weak 
coupling limit. The development of the method to the Bose gas is a subject of further 
works. 

With the notion 

:= e"^i ; D{t) := e^'^sD^-itHs 
the evolved interaction can be written in the form 

H,,,it) := Dit) ® + D+it) ® A+iS^g,)A{S^go) (1) 

The initial state of the compound system is supposed to be of the form 
P = Ps ® V^^^^- 

Here ps is arbitrary density matrix of the system and the initial state of the reservoir is 
the Gibbs state, at inverse temperature /?, of the free evolution, i.e. the gauge invariant 
quasi- free state ip^^\ characterized by 

= exp(-l </,(! + ee-^^^)(l - ee-^^^)-V >) (2) 

for each / G Tii. Here W{f) is the Weyl operator, j3 the inverse temperature of 
the reservoir, ^ = e^^ the fugacity, /i the chemical potential. We suppose that the 
temperature (3~^ > 0. Therefore for sufficiently low density one is above the transition 
temperature, and no condensate is present. The generalization to the case then the 
condensate is present, is a subject of further investigations. 

We will study the dynamics, generated by the Hamiltonian (|l]) and the initial state 
of the reservoir (|^) in the low density regime: n — > 0, t ~ 1/n (n is the density of 
particles of the reservoir). In the low density limit the fugacity ^ and the density of 
particles of the reservoir n have the same asymptotic, i.e. 

n^O n 

Therefore the limit n ^ is equivalent to the limit ^ — 0. 

The low density limit for this model, with completely different methods, based on 
quantum BBGKY hierarchy has been considered by Diimcke ||^. 

Throughout the paper, for simplicity, the following technical condition is assumed: 
the two test functions in the interaction Hamiltonian have disjoint supports in the energy 
representation. Thus the disjointness is invariant under the action of any function of H. 
This assumption means that the two test function go,gi in the interaction Hamiltonian 
satisfy: 

< go, S^e-^^gi >= Vt G M. 

This means that, even if the particles of the reservoir have generically a continuous 
energy spectrum, they behave like a 2-level system as far as their interaction with the 
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system is concerned: if Pq and Pi project onto disjoint intervals (energy bands) Iq and 
/i, these energy bands act as the counterpart of the energy levels eo, ei of the system. 

Also the rotating wave approximation condition will be assumed. This condition 
means that 

where uq is a real number. This is a familiar assumption, satisfied by all the 
Hamiltonians commonly used in quantum optics. This assumption is satisfied if 
D = |eo >< ei|, where |eo > and |ei > are eigenvectors of the free system Hamiltonian 
with eigenvalues cq and ei so that u = ei — eo. 

We will fix a projection operator Pq in Tii commuting with Hi and H and such 

that 

^05-0 = 9o and Pofl-i = 0. 
Using this projection let us define the group {St; t G M} of unitary operators on Hi by 

_ ^Og-iti^oPo _ ^it{Hi-uJoPo) 

The infinitesimal generator H[ of St is given by 
H[=Hi-uJoPo. 



Following Palmer |T3[ we realize the representation space as the tensor product 
of a Fock and anti-Fock representations. Then the expectation values with respect to 
the state ^p^^^ for the model with the interaction Hamiltonian (|I]) can be conveniently 
represented as the vacuum expectation values in the Fock-anti-Fock representation for 
the modified Hamiltonian. 

Denote by Hi the conjugate of Hi, i.e. 

6 : Hi^Hi, 6(A/) := A.(/) 

< ^(/),^(^) >,:= <g,f> 

then, H'l is a Hilbert space. The corresponding Fock space TlH'i) is called in this context 
the anti-Fock space. 

It was shown in that with notations Dq = D, Di = the modified Hamiltonian 
acting in T(Hi) ^{H'l) has the form 

H^{t) = 5^ ® {A+{Stg,)AiStgi-e) ® 1 

e=0,l 

+X{A{Stgi^e) ® A{StLg,) + A+{Stge) ® A+{StLgi^e)). 

Here A and A~^ are Bose annihilation and creation operators acting in the Fock spaces 
T{Hi) and 1(^0 and L := e'^^/^. 

Moreover, it will be assumed that there exists a subset IC (which includes go,gi) of 
the one particle Hilbert space Hi, such that 



/ \< f,Stg> \dt < oo Wf,gelC. 



Quantum stochastic equation 6 

The interaction Hamiltonian determines the evolution operator U^^^ which is the 
solution of the Schrodinger equation in interaction representation: 

with initial condition 

ut' = 1. 

One has the following integral equation for the evolution operator. 

t 



2. Energy representation 

We will investigate the limit of the evolution operator when ^ +0 after the time 
rescaling t t/^, where ^ = A^. After this time rescaling the equation for the evolution 
operator becomes 

dtUi^l = -i E ® (^^e,l-e,A(i) + i?l-.,.,A(^) + B+,_^^,{t))Ui^l 

£=0,1 

where we introduced the notations: 

Nsr,e,At) = ^A+{St/x^g,,)A{St/x2g,J ® 1 

Let us introduce the energy representation for the creation and annihilation 

operators by the formulae 

A+{g) = A+{P^g) = J dk{PEg){k)a+{k) = J dk5{H[ - E)g{k)a+{k) (3) 
AE{g) = A{PEg) 

Here 

oo 



-oo 



It has the properties 

St = j dEPE^'^ 

PePe' = S(E - E')Pe 
Pe = Pe 

In the case when Tii = L^(R'^) the one-particle Hamiltonian is the multiplication 
operator to the function u;{k) and acts on any element / e L^(R'^) as Hif(k) — 
oj{k)f{k). 
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It is easy to check that 

[AEif), A^g)] = S{E - E') < /, PEg > . 

Here < ■, ■ > means the scalar product in Hi. 
Using the energy representation one gets 

Be„e,At)= j dEidE2Be,,e,AEl,E2,t) 

where 

Qit{Ei-E2)/X'^ 

N,,^,,AEi,E2,t) := A+^{g,JAEM2) ® 1 

Qit{E2-Ei)/X-^ 

B,,^,,AEi,E2,t) := AE,{ge,)®AE,{Lg,^). 

Let us also denote 



je{E):= [ dt<g„Stge>e-''''. 



3. The limiting commutation relations 

Besides the operators i?£-^,£2_A(-Ei, -^2, ^) and A^ej,£2,A(-E'i, -E'2, ^) defined above let us 
consider the following operators: 

^U{E2-Ei)/\'^ 

K,e2AEi, E2,t) = l^A+^{Lg,,)AE2{Lgs2) 

with A~^{Lg^) has been defined in (^). The commutators of these operators are: 



[5ei,e2,A(-El, E2, t), 5+^^^ ^^(^3, E4, t')] 

Qi{t'-t)(Ei-E2)/X^ 

= [^eueA2,SiS{Ei - E3)6{E2 - E4) < ^^i, ^Ei^ei >< ge2, ^E2^^^e2 > 

+ X'6,,,,,5iEi - E,) < ge„PE,ge, > K,,e2AE,, E2, t') 

+ \^5e2,eME2 - E,) < g,„ PE2L'gs2 > iVs3,ei,A(^3, ^1, t')) (4) 

[-Bei,£2,A(-E'l, E2, t), A^£3,e4,A(-£^3, E4, t')] 

„i{t'-t)(Ei-E2)/X'^ 
= Ss,,s, ^2 ^(^1-^3) 

X < g,^,PE,ge, > fi,4,e2,A(^4,^2,t') (5) 
[Ne„e2AEu E2, t),N,,,,,AE3, E,, t')] 
e^' 



,i(t'-t)iE3-Ei)/X^ 



{5e2,eAE2 - E3) < g,^,PE2ge2 > Ne„s,AEi, E^, t') 



A2 

- 6e„eAEl - ^4) < ge„PE,gs, > iV.3,e2,A(^3, ^2, t) ) . (6) 
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Notice that in the sense of distributions one has the hmit 

iit'-t){Ei-E2)/\'^ 

lim = 2n5it' - - E2) (7) 

and, in the sense of distributions over the standard simplex (cf. 0) one has the hmit 

hm^ = 6+(t'-t)- = -. (8) 

A-.0 A2 +^ ^i(Ei -E2 -iO) ^ ' 

The following theorem describes the algebra of commutation relations for the master 

field in the LDL. 

Theorem 1 The limits 

X,,,,,{Ei,E2,t) := limX,,,,„A(^i,^2,t) {X = B,N) 

A — >0 

exist in the sense of convergence of correlators and satisfy the (causal) commutation 
relations 

[B,,^s,{E,,E2,t),Bl^^{E,,E,,t')] 

= 5e„eA2,eMt' - - E;)5{E2 - E,) 

^< 9ei, PEi9ei > ^ n t2 (n,\ 

^pr^£— io) 

[Bex,e2{El, E2, t), A^£3,j4(i?3, E4, t')] 

= 6,,Mt' - - ^3)||^z^^i?..,..(^4, E2, t') (10) 
[iV,,,,,(Ei, t), iVe3,.4(^3, E,, t')] = 6+{t' - t)— \ — 

H-C^3 — -fci — lUj 
X {5e2,eAE2 - E3) < gs2, ^^£2^.2 > N,^,eAEl, ^4, ^0 

-5,,,,,5{El - ^4) < 9er,PE^ge^ > N^.^.^iE^, E^, t)^ (11) 

The commutation relations of the master field are obtained by j^), i^TSj), (ff7|j replacing 
the factor 6+{t' - t) by 6{t' - t) and (i{Ei - E2 - 10))"^ by 27t6{Ei - E2). 

Proof. The proof of the theorem follows by induction from the commutation 
relations (^)-(§|) using the limits (0) and (^ and standard methods of the stochastic 
limit. 

4. The master space and the associated white noise 

Let Tio^i denote the closed subspace of Tii = L^(R'^) spanned by the vectors 

Stg„ e e {0, 1}, t e M. 
Let be a non zero subspace of Hq^ such that g^ & K {e = 0, 1) and 

< f, Stg > \dt < 00 \ff,geK. 
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This assumption implies that the sesquihnear form {-l-) : K x K — > C defined by 

oo 

{f\9)= j <f,Stg>dt, f,geK 

— oo 

is well defined. Moreover it defines a pre-scalar product on K. We denote {K, (-I-)} or 
simply K, the completion of the quotient of K by the zero (•|-)-norm elements. 

Define then Hilbert space Kq i which will be denoted also as K ®/3 K as the 
completion of the algebraic tensor product K ® K* with respect to the scalar product 

oo 

(/o®/5/i|/^®/3/0 := J <h,Stil><h,StL^f[><lt 

— oo 
oo 

= j <fo,Stfl>><f[,S.tL'h>dt. 



Bounded operators acts naturally on Kq^i by 

{A (^fs B){fo % /i) = Afo ®fs Bh VA, B G B{K). 

The limit reservoir (or master) space is the space 
T{L\^)®K^,i). 

The (non-causal) commutation relations (PD,..., (pJ]) mean that operators 



i?e;^,e2(-^i) -^2, ^) are the white noise operators ht{-) in JF(L (R) ® Kq 
The number operator is 

Ne,,,,iE^,E2,t) = J2 ne'{Ei)btige, 0/3 PE,ge'MPE,ge, ®pPE,ge') 

£' = 0,1 



where we denoted 



[ dEnAE,)BX,AE,E,,t)B,,,AE2,E,,t) (12) 

'=0,1 



< ge,PEL^ge >' 

By identifying the element of the algebraic tensor product f ^ g E K ^ K* with 
the operator 

\f><g\:^eni ^<g,^>feni 

so that 

(l/><^7|)* = l^7></l 
and introducing the scalar product of such operators X, Y (notice that = e~^^) by 

< F,X >:= Tr j e'^^Y* StX S;&t = 27rTr J dEe-I^^Y* PeXPe 
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one can rewrite white noise bt{g<^j3f) as bt{\g >< f\) with commutation relations defined 

by 

[bt{Y),b+{X)]^S{t'-t)<Y,X>. 
Let us introduce for simphcity 

Ne.AE.t) := j dE'N,^,,,iE,E',t) 
with (causal) commutation relations: 

[B,,,,,{E,t),Bl^^^{E\t')]^S,,,,,5,,Mt'-mE-E'^^^^^ > 
[B,,,, {E,t), iV,3,,, {E',t')] = 4„.3^+ {f - t) \e>-E-\^) ^^''^^ ^'^ • 

In these notations the hmiting Hamiltonian acts on Tis ® J^(-L^(R) ® -f^o,i) as 
H{t) = fdEj^De® {Ne,i-e{E, t) + B,.,,,{E, t) + Bt,,_,iE, t)). 

£=0,1 

5. Emergence of the drift term and annihilation process 

The results of the preceding section allow us to write the equation for the evolution 
operator in the stochastic hmit 

m = -iH{t)Ut = -iY,De® (iV,,i_,(t) + + BX,_M)Ut (13) 

£=0,1 

In order to bring it to the normally ordered form one needs to compute the commutator 



-iD,[Si_,,,(t),[/i] = -iD,j dE[B,_,,,{E,t),Ut]. 



Notice that D^Di^^ is a positive self-adjoint operator. Therefore one can assume 
that for each £■ e R, the inverse operator 

niE) (1 + {J,J^_,){E)D,D,_,)-' 

always exists. Notice also that, since Dg^Di^^ commutes with T^{E), one has 

1 - D,Di_,(7,7i_,)(^)r,(^) = T,{E). 

Therefore 

iL',(l - Di^,{j,ji^,){E)T,{E)D,) = iT,{E)D, (14) 

Theorem 2 For the model described above one has 
-iL>,[Si_,,,(i), = - J dE-f^_,{E)D,Dr_ME) 

X (< g,, PEL^ge >Ut- iD,^e{E)UtBi_,,,{E, t) + UtB,,,{E, t)) (15) 
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Proof. Using the integral equation for the evolution operator and the commutation 
relations (PD , (p!0|) , one gets 

-iD,[B^_,^,{E,t),Ut] 

t 

= -Y, DeDs' I dE' I dti[fii_,,,(E,t),iV,,i_,,(E',ti) + 5+i_,,(E',ti)]f/i, 

e'=0,l 

= -D,Z}i_,7i_,(E)(< ge^PEL^Qs > +B,^,{E,t))Ut (16) 

Notice that the first equality in (p!6| ) holds because, due to the time consecutive principle 

[5,,,,(E,t),f/iJ =0. 

Similarly one computes the commutator 

t 

[B,^,{E,t),Ut] = -iJ2De' J dE' J dt,[B,^,{E,t),N,,,,^,,{E' M)]Ut, 
^' 

= -iD,^,{E)B,^,^,{E,t)Ut (17) 
After substitution of this commutator into (|16D one gets 
-iD,[5i„,,,(E,t),f/i] = -D,Z}i_,7i_,(E)(< ge,PEL^ge > Ut-iD,^e{E) 

x([i?i_,,,(E,t), Ut] + UtB^_,,,{E,t)) + UtB,^,{E,t)) 

Then for 

UE,t):=-iD,[B^_,,,{E,t),Ut] 

one has 

(1 + (7,7i„,)(E)D,Di_,)/,(E,t) = -7i_,(E)D,Di_,(< ge,PEL^ge > Ut 

-iDME)UtBi-UE.t) + UtB,,,{E,t)). 

Since the inverse operator (1 + {^eli-e){.E)Di.Di^ey^ exists we can solve the 
equation above for fs{E,t). Using this solution we find 

-iD,[i?i_,,,(t),[/i] = j dEUE,t) = - j dE^^_,{E)D,D^_,T,{E) 

X (< gs,PEL'gs >Ut- iD,^,{E)UtB,_,,,{E,t) + t/,5,,,(E,t)). 

6. Emergence of the number and creation processes 

In order to bring equation ([13|) to the normally ordered form one needs also to move 
the annihilation operators in A^^ i_£(t) to the right of the evolution operator. Using (0) 
this leads to 

- iD,iV,,i_,(t)t/i = -iD, [dEj2 n,,{E)B^^^,{E,t)B^_,^,,{E,t)Ut 

e'=0,l 

= -iD, I dEY, n,>{E)B^^^,{E,t)([Bi^,,,iE,t),Ut] + UtBi.e,e'{E,t) 

^ e'=0,l 

and the commutator is evaluated using (|16D and (|l^). 
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Theorem 3 For the model described above one has 

+ nAE){iT,iE)D,B^^^,iE,t)UtB,.,AE,t) 

e'=0,l 

+-fi.e{E)D,Di^,T,{E)B+^,{E,t)UtB,,,>{E,t)')^ 
Proof. From (|16]) and ([T7|) it follows that 

-iD,^,iE)UtB^^,,AE,t) + UtB,,AE,t)). 



After substitution of these commutators in (|T8|) and using ( [I^ ) one finishes the proof of 
the theorem. 

7. The normally ordered equation 

Theorems (0) and allow us to obtain immediately the normally ordered equation for 
the evolution operator in the LDL. This procedure of deduction of quantum stochastic 
differential equation is being called a stochastic golden rule. The normally ordered 
equation has the form 

dtUt = -J2f dE{iT,{E)D, J2 nAE)B+^,{E,t)UtB,^,AE,t) 

£=0,1 £'=0,1 

+7i_,(i?)D,Di_,T,(E) J2 nAE)B+^,{E,t)UtB,^AE,t) 

£' = 0,1 

+iT,{E)D,B+,_^{E,t)Ut + -fi-e{E)D,Di^,T,{E)B+^{E,t)Ut 
+iT,{E)D,UtB,^,,,{E,t)+j,^,{E)D,D^_ME)UtB,,,{E,t) 
+-i,.,{E)D,D,_,T,{E) < g,, PeL^Qs > Ut} (19) 



We will represent equation (|I9|) in the form of quantum stochastic differential 
equation |]T^ . 



Theorem 4 Equation is equivalent to the quantum stochastic differential equation 
dUt = -J2 I d^{ir,(E)I),diV,(27r|(7, >< g^_,\PE ®p Pe) 

. n 1 "J 



£ = 0,1 ' 



+^^.,{E)D,D^_,T,{E)dNt{2'K\g, >< g,\PE ®/3 Pe) 
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+iT,{E)D,dB+{g, ®p PEgi-e) + ^i-e{E)D,Di_,T,{E)dBt{ge ®p PEQe) 



+iT,{E)DABt[gi-e ®p Pe9s) + ii-e{E)D,D^.,T,{E)dBt{ge ®p PEQe) 
+^,.,{E)D,D,_,T,{E) < g,, PEL'g, > dt^Ut (20) 

Proof. Let us consider the following term in (|T^ 

J2 n,,{E)B+^,{E,t)UtB^.,AE,t) = n,_,{E)bt{ge ®/3 PE9i~e)Utht{gi^e ®p PeQi-s) 

£' = 0,1 

+n,{E)ht{ge ®p PEge)Utbt{gi-e ®/3 PeOs) (21) 

The matrix element of this expression on the exponential vectors ip{f), ip{f') is (we use 
Dirac's notations also for bra- and ket-vectors from Kq^i) 

(ni„,(E) < f\g, PeQi-s >< 9i-e ®p PEQl-elf > 

+n,iE) < f\g, ®^ PEge >< gi-e ®P PE9e\f' >) < ^1)^^/') > 

=< /|T,,i.,(E)i/' X ^if)\utmn > (22) 

which is the time derivative of the matrix element < il!{f)\dNt{T^^i_^{E))\ilj{f') > where 
Nt{Ts^i-s{E)) is the number process with intensity 

Te,l-e{,E) = ni^,{E)\gs ®p PeQi^s >< 9l-e ®/3 PeOi-sI 
+neiE)\ge 0/3 PEge >< gi-e ®/3 PeQe] 

Let us now prove that 

T,^i-e{E) = 27i\g, >< g,_,\PE ®/3 Pe (23) 
To this goal let us consider the action of the Ts^i-s{E) on vectors of the form 

1/ >= \PE,9ei ®pPE29e2 > ■ 

One has 

T,,i_,(E)|/> 

= 27TSe,l-eAEl - E)S{E2 - E) < gi^„ PEgi-e > 
(4,l-£2lfi'£ ®/3 ^i?5'l-£ > +4,£2lfi'£ ®/3 PE9e >) 

= 27r5,,i_,, < (7i-£, PE9i-e > KEi - E)6iE2 - E) 

{5e,l-e.2\9e ®/3 ^Efe > +4,£2lfi'£ ®/3 ^Sfe >) 

= 27r5,,i_,, < PEgi-e > KEi - E)6{E2 - E)\g, Psgs^ > 
= 2Tx\ge><gi-e\PE®pPE\f > (24) 
Therefore (^) holds and the term ( plj) corresponds to the number process 

dNt{2Tx\ge >< gi-e\PE ®/3 Pe) 
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Computing the same matrix element for the term 
J2 nAE)B+^,{E,t)UtB,^AE,t) = n,_,{E)bt{g, ®^ PE9i-e)UA{g, ®p PE9i-e) 

e'=0,l 

+n,{E)b+{g, ®^ PEge)Utbt{gs ®/j P^^e) (25) 
one finds an expression hke (22) with i_£ replaced by the operator 

Te,e{E) = ni_e{E)\ge ®p ^E^l-e >< ®/3 ^E^l-el + ne{E)\ge ®p PEge >< 9e ®/3 P_B^e| 

and a calculation similar to the one done in (^^ leads to the conclusion that 

T,^,{E) = 27r\g,><g,\PE®^PE. 

Therefore the term (^) corresponds to the number process 

dNt{27T\g, >< ge\PE®^PE). 

This finishes the proof of the theorem. 
Let us introduce the notations: 

R,,e{E) = --fi-,{E)D,Di^,T,{E) 
Rs,i-e{E) = -iT,{E)D,. 

In these notations the quantum stochastic differential equation for the evolution operator 
can be rewritten as 

dUt=J2 I dE^Re,i-e{E)dNt{2'K\g, >< g^_,\PE ®p Pe) 

e=0,l 

+Re,e{E)dNt{2Tl\ge ><ge\PE®p Pe) 

+R,,^.,{E)dBt{ge ®/3 PeQI-s) + ReAE)dBt{ge ®/3 PEge) 
+Re,l-e{E)dBt{gi-e ®/3 PEge) + Re,e{E)dBt{ge ®/3 PEge) 

+Re,e{E) <gs,PEL^ge>dt}Ut. 

Notice that the quantum stochastic differential equation ( pOD can be written also in 
the Frigerio-Maassen form [0. In order to prove this recall that for any pair of Hilbert 



spaces Xq, Xx if A^, A denote the number and annihilation processes on the Fock space 



J^{Xi) then for Xq G B{Xo), Xi G B{Xi), x G Xi, Frigerio and Maassen introduced 
the notation: 

N{Xo®Xx) :=Xo®iV(Xi) 

A{Xo ® Xix) := Xq ® A{Xix) 

< x,Xo® Xix >:= Xo®l < X, Xix > 
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Let us also introduce an operator T^{E) acting on the triple 71$ ® K ®p K (this is the 
reason for introducing index 3) as 

T^{E):=2lT J2 ReAE)®\9e><9e'\PE®pPE 



e,e'=0,l 

and the vector ^{E) e K ®p K 



i{E) := 7^ ^ -9e ®p 9e- 

< 9e, Pe9s > 



In these notations equation (^Dj) can be written as 
dUt = J dE(dNt{n{E)) + dB+{n{E)^{E)) 

+dBt{T;{E)aE))+ < aE),T3{E)C{E) > dt^t- 

8. Connection with scattering theory 

Here we consider relation between the evolution operator and scattering theory. Because 
of the number conservation, the closed subspace of TCs ® J-^ generated by vectors of the 
form {u e Hs, f e Hi = L'^{W^)) which is naturally isomorphic to Hs®Hi, 

is globally invariant under the time evolution operator exp[i(ifs (8> 1 + 1 (8> -f^R + V)t]. 
Therefore the restriction of the time evolution operator to this subspace corresponds to 
an evolution operator on ® Hi given by 

exp[i(ifs ®l + l®Hi + Vi)t] 

where 

Vi= J2^-®\9e><9i-e\ (26) 

£ = 0,1 

The 1-particle M0ller wave operators are defined as 
n± = s- lim exp[i(i^s ®1 + 1® H + Vi)t] exp[-i(i^s ® 1 + 1 ® H)t] 
and the 1-particle T-operator is defined as 

T = 1^1^+ (27) 
From (p6|) it follows that 

n± = s- lim Ui'^'^ 

t— >itoo 



where Uj; is the solution of 



(1) 

dtUi'^ = -i ( 5^ ® \St9e X St9l-e\ ) Ui'^ f/^') = 1. 

Ve=0,l / 

In order to make a connection between the stochastic process Ut and scattering 
theory notice that the operator T^iE) can be written as 

Ts{E) = 27iT{E) ®p Pe 
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where operator T{E) acts on Hs ® K as 

T{E)= J2 Re,e'{E)®\g,><g,,\PE. 

£,£'=0,1 

In it was proved that T-operator defined by ( p?] ) connected with T{E) by the following 
formula 

dET{E). 



9. Vacuum expectation value 

For the vacuum matrix element of the evolution operator from (^) one immediately 

gets 

< U{t) >,ac= e-^*. (28) 
The operator T acts in Hs as 

r = E / dEj^_,{E)D,D^_,T,{E) < ge,PEL'ge > ■ 

£=0,1 

Theorem 5 The operator T has a non-negative real part (i.e. this operator describes 
the damping). 

Proof. From the definition of T^{E) we know that 



^,^,{E)D,D^.,T,{E) 



7r-£(i?)+7£(^)^£A-£ 
But 7e and 7i_£ (hence also l/7i_e) have positive real part and DeDi^s is positive 
self-adjoint. 

Hence the above expression has a positive real part because it is of the form: 
H _ H(Re zi + H Re Z2) .iy(lm Zi + Hlui Z2) 



Z1+Z2H \Z, + Z2H\^ \ZI+Z2H\^ 

where H is positive self-adjoint and Zi, Z2 have a positive real part. Since PEL'^ge) > 
the thesis follows. 
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